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Real-time thermal field theory is known in the two flavors “closed-time path formalism” and 
“thermo field dynamics”. With a certain choice of parameters the full two-point functions of these 
models are identical, hence a scheme to diagonalize the Green function in one model can be applied 
to the other as well. This allows to compare the diagonalization schemes that have been discussed 
in the recent literature in order to select the simplest one. Such comparison speaks in favor of a 
diagonalization method which is motivated from thermo held dynamics. 


In real-time quantum field theory for statistical systems, i.e., for thermal as well as for non-equilibrium states, the 
two-point Green functions acquire a 2 x 2 matrix structure. Two flavors exist for such a theory: The Schwinger-Keldysh 
closed-time path formalism (CTP, [Q) and thermo field dynamics (TFD, 0). 

These matrix valued Green functions (propagators) contain spurious information that unnecessarily complicates 
their use in a perturbative expansion. This unneeded information may be removed by matrix diagonalization. The 
fact that different diagonalization methods are possible has led to a dispute in the recent literature ; the present 
work is aimed at a clarification of this question, thereby uniting the different approaches. 

To this end, we concentrate on an interacting scalar quantum field in a thermal equilibrium state and note that 
similar considerations apply for fermions. The first question to be settled is, whether the Green functions under 
dispute are different. By definition, the propagator matrix in the GTP formalism is |Q] 



(T [4>a;(t>x'] ) {(jix'fpx) \ 

{(j^xfl^x') (T[4 ’x4‘x>]^ j 


( 1 ) 


where (■) denotes the statistical average, T [•] the time ordered and T [■] the ’’anti-time ordered” product, i.e., 

T [4ix4'x>] = 0(a;° - x°') ipx'px' + - x°) (j)x>(l>x 

T [(j^xfkx'] = 0(a;°' - a;°) + 0(a;° - x°') . (2) 

In the TFD formalism, the canonical commutation relations have two different commuting representations and pxi 
the matrix valued propagator is 




(T [pxfpx'] ) 




( 3 ) 


TFD is conceptually different from GTP, in the sense that it contains a Bogoliubov transformation of quantum fields 
at the operator level. This transformation contains three parameters, for the purpose of the present paper only one 
aspect is interesting: For a certain choice of TFD parameters, the two matrices above are identical |^|. By adopting 
this special choice (a = I), we will henceforth omit a distinction between Dqtp and D^pp), and conclude that they 
may be diagonalized by the same matrix transformation. 

We now turn to the task of isolating the spurious information hidden in these propagators. For convenience, 
we perform a Fourier transformation to the momentum variable p = {po,p), with respect to the difference of the 
space-time coordinates. 

The goal is to formulate a perturbative expansion in terms of the retarded and advanced functions 

D^{p) = D^^ip) - D^^p) 

D^{p) = D^\p) - D^\p) = {D^ip)y , (4) 
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since they are free of poles on the physical complex energy sheet other than on the real axis. By inspection of (|^) and 
one may easily realize that 

D^^p) + D^^ip) = D^^p) + D^^p) 

D^^p) =-{D^\p)y 

D^\p) = -{D^yp)y , D^\p) = - {D^\p)y , ( 5 ) 

where * denotes complex conjugation. Apart from these trivial relations we also know, that in a thermal equilibrium 
state the Kubo-Martin-Schwinger boundary condition holds, 

{l + nB{po))D^'^(p) - nB{po)D'^^{p) = 0 , ( 6 ) 

with the Bose-Einstein function 

= exp(/3(£; -p))-l 

at inverse temperature (3 and chemical potential p {p = 0 for the real scalar held). 

From these relations follows that there exists a nonsingular 2x2 matrix V{pq) such that 

Vipo) Dip) Vipo) = ( ) , (8) 

where 

D^p) = 0(po) D^ip) + 0(-po) D^ip) . (9) 

(this “Feynman”-like propagator is not causal in a state with nonzero nB(po))- 

This matrix transformation diagonalizes the full boson propagator. The transformation matrix is well known for 
some time (see e.g. [Q, eq. (2.4.31)]), here we adopt the notation of ref. [§: 

Vfr, \ - ( cosh( 6 l) - exp[/3(po - m)/ 2 ] sinh( 6 »)) \ , . 

(Po) \^ — exp[—/3(po — At)/ 2 ] sinh( 0 ) cosli( 0 ) J ’ ^ ' 

with hyperbolic functions 

cosh( 0 ) = \/ 0 (po)[l+nB(po)] - 0 (-Po)?^B(po) 

smh( 6 ») = v^ 0 (po)nB(po) - 0 (-Po)[l + ?^b(po)] • ( 11 ) 

Thus, after diagonalization the only temperature dependence is located in the matrices V - and we may absorb 
them into the interaction vertices of a theory, which are then linked by retarded or advanced propagators. The goal 
of isolating the spurious information hidden in the propagator from its field theoretical content has been reached. 
However, three questions are unanswered yet: 

1. Does the matrix represent the simplest possible choice ? 

2. What is the physical meaning of this diagonalization procedure ? 

3. How does one handle the diagonalization task in non-equilibrium states ? 

To answer the first question, we reconsider the relations (^ and (^). Apart from the diagonalization procedure 
outlined above, they also guarantee a different way to reach our goal: There exists at least one nonsingular 2x2 
matrix B such that 

B[nBiPo)] Dip) T3 ^H[nB(po)]^ ^ ^ D^ip) ) ’ 

where T 3 = diag(l, —1). Hence, also this matrix transformation diagonalizes the full boson propagator 

In fact, more than one such matrix B exists, but one may chose a particularly simple form linear in the Bose-Einstein 
function n(po) 
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„r , , (1 + n) n 

= 1-1 1 


( 13 ) 


From the continuum of possible diagonalization methods of the full propagator, the above choice clearly is one of the 
simplest - whereas the matrix V in (0) is a highly nonlinear function of n(po). 

To answer the second question, we consider a simple harmonic oscillator with hamiltonian iJ = u; a, immersed 
in a heat bath of inverse temperature /3. It is well known, that one may describe this oscillator equally well in terms 
of particle states or hole states, i.e., the Liouville operator governing the time evolution of the statistical system has 
a symplectic symmetry. Without elaboration at this point we note that its symmetry group is the two-dimensional 
symplectic group Sp(2) (see 0i for a more complete discussion). 

In a thermal equilibrium state, we also know the statistical operator to be VF = exp(—/3FI). It follows that 


(1 + nB(ijj)) aW —nB(i^)Wa = 0 
(1 -I- nB{(^)) W — nB(i^) W= 0 ^ 


(14) 


which is nothing but the Kubo-Martin-Schwinger condition for this oscillator. Thus, although the creation and 
annihilation operators of this simple system do not commute with the statistical operator, we find that a certain 
linear combination of right-acting and left-acting operators (well defined in Liouville space [p^ ) annihilates the 
density matrix. 

The corresponding orthogonal linear combination creates an excitation of the system. It is a linear combination 
of particle and hole state, which propagates through the system without “feeling” the thermal background. Thus, 
by introducing these “thermal quasi-particles” the statistical information has been separated from the time evolution 
problem. 

The relation between ordinar^reation/annihilation operators and these unusual linear combinations is mediated by 
the matrices B\n{u!)\ from eq. (p^). Hence we may identify the diagonalization according to eq. (0) physically as the 
transformation to a linear combination of particle and hole state that have retarded/advanced boundary conditions 
in time. 

It remains to answer the third question, i.e., what one does for non-equilibrium states. In this case, the Fourier 
transform with respect to coordinate differences leaves us with an additional dependence on X = {x + x')/2. In this 
mixed representation, the relations i) between the matrix elements of the two-point functions still prevail, but the 
KMS boundary condition d) is not valid. 

We label the non-equilibrium propagators by a~-sign, i.e. 1?“*' = D°'^{X,p), and use the easily established fact that 
for arbitrary parameter N = N{X,p) the matrix B[N] according to 0) leads to 


B[N]D{X,p)t3 




7V)D21 

]jR 


(15) 


Hence, this procedure transforms the full Green function to triangular form. One might then ask, whether there exists 
a special choice for N = N{X,p) such that the off-diagonal element of this matrix also vanishes, i.e. such that 


N{X,p)D^\X,p) - (1 + NiX,p)) D^^X.p) = 0 . 


(16) 


This can be answered by inserting the propagator into the full Schwinger-Dyson equation, the result is a differential 
equation for N(X,p). On close inspection it is identified with a transport equation, i.e., a close relative of the 
Boltzmann equation for the quantity N{X, p) [p|JTI|] . 

Naturally, transport equations are also obtainable using a different diagonalization scheme. However, the differential 
equations then are much more complicated and thus their relation to the Boltzmann equation or any other known 
transport equation is not easily established. 

To conclude the present paper: As already stated in ref. |^, a diagonalization of two-point functions is easily possible 
in the closed-time path formalism as well as in thermo field dynamics - which is a trivial fact since the propagators 
can be made identical. 

As pointed out in ref. as well as in Q, the diagonalization matrices that occur naturally are V from (0) in CTP, 
but B from 0) in TFD. For simplicity one should therefore make the obvious choice of the diagonalization scheme 
according to eq. ( |l^ also in the CTP formalism. Especially for non-equilibrium states such simplicity is required to 
obtain a meaningful interpretation of the diagonalization condition - independently of the reader’s preference for one 
of the two flavors of real-time statistical quantum field theory. 

The situation described here bears a close analogy to the fixing of a gauge: In equilibrium states, the choice 
of diagonalization matrices does not influence physical quantities - but the calculational effort greatly depends on 
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this choice. It has been shown, that this is indeed more than an analogy, i.e., that choosing a special form of the 
transformation matrices indeed corresponds to a gauge fixing ^ . 
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